A Systematic Study of the Superconducting and Normal State Properties of Neutron 

Irradiated MgB2 



O 

o 



(N 
(N 

o 
o 

u 

Oh' 



R. H. T. Wilke, S. L. Bud'ko, and P. C. Canfield 

Ames Laboratory US DOE and Department of Physics and Astronomy, Iowa State University, Ames, lA 50011 



J. Farmer 

Missouri University Research Reactor, University of Missouri 



Columbia, Research Park, Columbia, MO 65211 



S. T. Hannahs 

National High Magnetic Field Laboratory, Florida State University, 
1800 E. Paul Dime Drive, Tallahassee, Florida 32310 
(Dated: February 2, 2008) 

We have performed a systematic study of the evolution of the superconducting and normal state 
properties of neutron irradiated MgB2 wire segments as a function of fluence and post exposure 
annealing temperature and time. All fluences used suppressed the transition temperature, Tc, 
below 5 K and expanded the unit cell. For each annealing temperature Tc recovers with annealing 
time and the upper critical field, Hc2(T=0), approximately scales with Tc. By judicious choice 
of fluence, annealing temperature and time, the Tc of damaged MgB2 can be tuned to virtually 
any value between 5 and 39 K. For higher annealing temperatures and longer annealing times the 
recovery of Tc tends to coincide with a decrease in the normal state resistivity and a systematic 
recovery of the lattice parameters. 

PACS numbers: 74.25.Bt; 74.25Fy; 74.25.Ha 
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I. INTRODUCTION 

It has been well established that superconductivity 
near 40 Ki in MgB2 is phonon mediated^*^ with two dis- 
tinct energy gaps'* [see reference 5 and works therein for 
a review of the basic properties of MgB2]. Much recent 
work has focused on understanding how structural de- 
fects and chemical substitutions affect the superconduct- 
ing properties. The goal is two fold; first, to understand 
the basic physics in this unique system and second, to 
use this understanding as a basis for developing practical 
superconducting wires, tapes, or devices. 

Of the many elements investigated, only Al and C have 
shown concrete evidence of entering the structure^. Both 
electron dope the system and it is the the Fermi sur- 
face changes, rather than interband scattering, which are 
believed to cause the suppression of T^—il. When alu- 
minum substitutes for magnesium, Hjlj'' decreases while 

■p"_Lab 
J^c2 



remains constant or slightly increases, resulting in 



an decrease in the anisotropy ratio, ^ h—^\2' 
These changes can be understood as the result of changes 
in the Fermi surface topology of the a bandit. In the case 
of carbon substitution for boron, the upper critical field 
is enhanced in both directionsii*i^, with H;f^°'' increas- 
ing more rapidly than Il|!2'' leading to a decrease in the 
anisotropy ratioi. This enhancement of the upper criti- 
cal field can not be explained in terms of changes in the 
Fermi surface and is believed to result from an increase 
the scattering in the tt band^ in accordance with the the- 
ory proposed by Gurevichi^. 

Damaging or disordering of the sample using protons, 
heavy ions, neutrons, etc., is another route to systemati- 
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cally changing the system. Of these possible routes, neu- 
tron irradiation offers the best avenue for uniformly dam- 
aging bulk MgB2. There are two main sources of damage 
from neutron irradiation of MgB2. First, fast neutrons 
deposit energy through inelastic collisions with atoms, 
creating thermal and dislocation spikesi^. Second, *°B 
has a large capture cross section for lower energy neu- 
trons and readily absorbs these thermal neutrons, subse- 
quently a decaying to ^Li. Early neutron damage stud- 
ies focused on irradiation of powders or pressed pellets of 
MgB2 containing natural boron-'^^*iSiiLi&. 

The absorption of slow neutrons by ^"B can lead to 
self shielding and prevents uniform damage in bulk MgB2 
samples significantly larger than the penetration depth. 
Under the assumption of linear absorption, the intensity 
of the incident beam decreases exponentially within the 
sample, and is given bj*i^: 



/(<) = Ioexp{-tsN) 



(1) 



where Iq is the initial intensity, t is the sample depth, 
s is the absorption cross section, and N is the number 
density of atoms. Here, N corresponds to the number 
of *°B atoms per cm'^. A calculation of the half depth 
in MgB2 synthesized from natural boron, which consists 
of 19.9% ^"B, yields a depth of approximately 130 /im. 
Samples with dimensions significantly larger than 130 
/im will contain substantial gradients of damage associ- 
ated with slow neutron absorption. Therefore several dif- 
ferent approaches have been employed to minimize neu- 
tron capture, ensuring homogeneous damage throughout 
bulk samples. Thermal neutrons have been blocked by 
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a cadmium shiel(iiSi22i2i or the natural boron has been 
replaced with isotopically enriched ^^BSSi^S.. 

An alternate approach is to damage MgB2 objects with 
characteristic dimensions comparable to or smaller than 
130 /zm. In this case natural boron can be used and 
relatively homogeneous absorption of slow neutrons can 
be assumed. In this paper we report on the isotropic 
irradiation of fully dense, 140 /xm diameter, MgB2 wires 
with both thermal and fast neutrons. Although these 
fibers contain natural boron, they were exposed to an 
isotropic fluence and had sample dimensions comparable 
to the penetration length, implying the defect structure 
should be fairly uniform throvighout the samples. 



II. EXPERIMENTAL METHODS 

We synthesized fully dense, 140 micrometer diameter 
MgB2 wire by reacting boron filaments with magnesium 
vapor at 950 °C for 36 hours. This technique is described 
in full detail elsewhere^'*. Once synthesized, three fibers, 
each approximately 1-2 cm in length, were sealed under 
a partial helium atmosphere in quartz ampoules with di- 
mensions of 3 mm I.D., 4 mm O.D., and approximately 
3 cm long. A partial atmosphere of inert gas was neces- 
sary to provide thermal contact with the cooling water, 
so as to prevent the filaments from overheating during 
the irradiation. Helium was used due to its low neutron 
capture cross section. A 2.5 cm diameter, 4.7 cm long 
water flooded aluminum can containing 25 ampoules was 
then exposed to an isotropic flux of reactor neutrons, con- 
sisting of 98% thermal neutrons (E=25.3 meV) and 2% 
epithermal neutrons (ranging in energy up to 10 keV), 
at the Missouri University Research Reactor (MURR) 
for time periods of 24, 48, 72, and 96 hours. The expo- 
sures corresponded to fluences of 4.75 x 10^*, 9.50 x 10^^, 
1.43 X 10l^ and 1.90 x 10^^ cm^^ respectively (Table 1). 
Table 1 also includes estimates of the atomic percentage 
of boron converted to lithium. The fluences were esti- 
mated using a cobalt doped aluminum flux wire. The 
percent of lithium produced was estimated, under the 
assumption of a uniform fluence throughout the sample, 
by the product of the density of ^°B, the ^°B neutron 
capture cross section, and the neutron fluencaiS,. Post 
exposure anneals were performed with the samples still 
sealed within the quartz ampoules for temperatures up 
to 500 °C in a Lindberg model 55035 Mini-Mite tube fur- 
nace. In each case, we annealed an ampoule at a given 
temperature for a set period of time. Upon removal from 
the furnace, the ampoules were quenched in air and then 
opened in order to perform measurements on the individ- 
ual wires. Each data point (i.e. each fluence/anncaling 
time/annealing temperature) is from a single sample that 
was annealed only once. In some cases we measured sev- 
eral samples which underwent nominally the same con- 
ditions and this was done either by using different wires 
within a given ampoule or from annealing additional am- 
poules under the same profile. In all cases measurements 



were performed on samples which underwent only a sin- 
gle anneal. 

The boron filaments used in this study contain a tung- 
sten boride core. Fast neutrons colliding with ^^^W 
atoms, which have a natural abundance of 26.3%, can 
be absorbed into the nucleus causing the emission of a 
proton and transforming the tungsten into -'^^^Ta. ^*^Ta 
f3 decays back to ^^^W, with a half life of 181 days. As a 
result the filaments were mildly radioactive and required 
appropriate safety measures in handling. 

Powder x-ray diffraction (XRD) measurements were 
made at room temperature using CuKa radiation in a 
Rigaku Miniflex Diffractometer. A silicon standard was 
used to calibrate each pattern. The experimentally de- 
termined Si peak positions were found to be offset from 
their known values by a constant amount. Within each 
spectra, the peaks varied about some constant offset and 
this variation was used to estimate experimental uncer- 
tainty in the lattice parameters. Lattice parameters were 
determined from the position of the (002) and (110) 
peaks. DC magnetization measurements and magneti- 
zation hysteresis loops were performed in a Quantum 
Design MPMS-5 SQUID magnetometer. Transport mea- 
surements were done using a standard AC four probe 
technique, with platinum wires attached to the samples 
with Epotek H20E silver epoxy. In a previous study— 
we note that during the synthesis of the boron filament, 
the tungsten reacts with B in the gas stream forming a 
host of tungsten boride compounds. As a result the core 
is not a low resistance W wire, but a much more highly 
resistive alloy of the refractory boride materials W^^Bj, 
that, due to expansion associated with the formation of 
the MgB2 is segmented. In addition to the transmutation 
of some of the W to Ta, the presence of natural boron 
in the core will lead the presence of defects caused by 
the transmustion of ^°B to ^Li just as in the MgB2 wire. 
Additionally, the cross section of the core is only 2% of 
the cross section of the wire. We therefore take the mea- 
surement of the resistivity to represent the resistivity of 
the MgB2 . (This assumption is valid as long as the resis- 
tivity of the core remains above ~ O.bfifl cm. Given the 
non-stoichiometry, radiation damage, transmutation and 
segmentation of the core we feel this is a safe assump- 
tion.) Resistivity versus temperature measurements in 
applied magnetic fields up to 14 T were carried out in a 
Quantum Design PPMS-14 system and resistivity versus 
field was measured up to 32.5 T in a DC resistive magnet 
using a lock-in amplifier technique at the National High 
Magnetic Field Laboratory in Tallahassee, Florida. 



III. FLUENCE AND ANNEALING STUDIES 

Normal state and superconducting properties were 
found to be a function of fluence, annealing time, and 
annealing temperature. We present the data in the fol- 
lowing order: (a) structural studies on as-damaged wires, 
(b) 24 hour anneals at varying temperatures for a flu- 
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ence of 4.75 x 10^* cm~^, (c) Variable time anneals for 
a fluence of 4.75 x 10^^ cm"^, (d) 24 hour anneals for 
all measured fluence levels, and (e) 1000 hour anneals 
on all four fluence levels. With these cuts through the 
multi-parameter, fluence, annealing temperature, anneal- 
ing time phase space we can start to describe the effects 
of neutron damage on MgB2 samples made with natural 
boron. 



A. Structural Studies of As Damaged Samples 

The x-ray (002) and (110) peaks for an un-damaged 
wire and the entire set of as-damaged samples are plot- 
ted in figure la. The lowest exposure level, 4.75 x 10^^ 
cm~^, shows an anisotropic expansion of the unit cell. 
The a- lattice parameter increases from 3.0876(5) A in 
the un-damaged sample to 3.0989(2) A, an increase of 
0.0113(7) A or 0.37%. The c- lattice parameter increases 
from 3.5209(7) A to 3.5747(2) 1, an increase of 0.0538(9) 
A or 1.02%. Similar anisotropic expansion of the unit 
cell was seen by Karkin et. al»i^. The authors report lat- 
tice parameter increases of Aa=0.0075 A or 0.24% and 
Ac=0.0317 A or 0.9% for a fluence of 1 x 10^^ cm'^ 
thermal neutrons and 5 x 10^* cm~^ fast neutrons. For 
irradiation of isotopically enriched Mg^^B2, little change 
was seen in the a- lattice parameter up to a fluence level 
of 10^^ cm~^— . For this fluence level the authors report 
a 0.008 A or 0.23% increase in the c- lattice parameter 
relative to an undamaged sample. Due to the large dif- 
ferences in Aa and Ac reported for different irradiation 
conditions, we performed systematic study of the lattice 
parameters as a function of fluence level, annealing tem- 
perature, and annealing time for samples all prepared in 
the same manner. 

Relative to the 4.75 x 10^* cm~^ fluence level, filaments 
exposed to a fluence of 9.50 x 10^^ cm^^ show a fur- 
ther expansion of the unit cell, with the a- and c-lattice 
parameters reaching 3.1017(7) A and 3.5805(10) A, re- 
spectively. The response to further increases in exposure 
qualitatively changes. For an exposure level of 1.43 x 10^^ 
cm~^ the (002) peak broadens substantially and the a- 
and c- lattice parameters contract relative to the samples 
exposed to the 9.50 x 10^* cm~^ fluence level. Increasing 
the fluence further to 1.90 x 10^^ cm~^ results in a further 
contraction of both the a- and c- lattice parameters with 
the (002)peak remaining broad. The x-ray scan for the 
as-damaged, 1.90 x 10^® cm~^ fluence level sample, from 
26'=20-65° is plotted in flgure lb. The peak widths for 
the MgB2 (hkO) peaks are within a factor of two of the 
widths for neighboring Si peaks, whereas the MgB2 (001) 
peak is not resolvable and the (002) peak is more than 
flve times as broad as the Si peaks. It should be noted 
that in pure, undamaged MgB2 all x-ray peaks have full 
width at half maximum (FWHM) values comparable to 
those of the Si standard. For example, the peaks used to 
determine the lattice parameters in the undamaged wires 
used in this study, the (002) and (110), had FWHM val- 



ues of 0.216° and 0.221° respectively. The neighboring 
Si peak, the (311) peak which occurs at 56.110° 26, had 
a FWHM of 0.190°. 

In the case of a fluence level of 1.90 x 10^^ cm~^, 
the MgB2 (100) peak, occurring at 261=33.397° has a 
full width at half maximum (FWHM) of 0.216°, in com- 
parison to the FWHM of 0.180° for the Si (111) peak 
at 261=28.422°. Similarly the MgB2 (110) and Si (311), 
which occur at 59.686° and 56.110° 26*, had FWHM val- 
ues of 0.291° and 0.163° respectively. In contrast, the 
MgB2 (002) peak has a FWHM of 1.335°. The resul- 
tant correlation lengths, given by A29 = X/d cos{9), are 
approximately 4000 A in the plane and 700 A along the c- 
direction. The mixed (101) peak has a FWHM of 1.045° 
and a correlation length near 1000 A, values intermedi- 
ate between the two extrema. The two highest fluence 
levels have resulted in an anisotropic change in the corre- 
lation length, with a decrease in long range order between 
the boron planes. Since the (002) peak has substantially 
broadened for higher fluence levels, tracking the trend 
in Ac becomes harder due to the inherent decrease in 
accurately defining peak position. 



B. 24 Hour Anneals of Samples Exposed to a 
Fluence of 4.75 x 10^** cm"^ 

As-damaged samples showed a suppression of Tc to 
below 5 K. This is a much larger suppression than seen 
upon irradiating Mg^^BjS^: the authors report a Tc of 
12.2 K for a fluence of 3.9 x 10^^ cm~^. Irradiation of 
MgB2, i.e. with natural boron, using a fast neutron flu- 
ence of 2 X 10^8 cm-2 gave a T^ of 30.2 Kii. Thus dam- 
age induced via neutron capture and subsequent alpha 
decay appears to play a signiflcant role in suppressing 
superconductivity. To further investigate the effects of 
uniform irradiation on samples containing natural boron 
we have performed a systematic study of the normal state 
and superconducting properties as a function of anneal- 
ing temperature for samples annealed for 24 hours. 

We annealed a set of wires, exposed to a fluence of 
4.75 X lO^* cm-^ for 24 hours at 100 °C, 150 °C, 200 
°C, 300 °C, 400 °C, and 500 °C. X-ray measurements 
indicate the initial expansion of the unit cell could be 
systematically reversed by subsequent annealing (Figure 
2), with the Aa and Ac values decreasing with increas- 
ing annealing temperature (Figure 2b). The a-lattice pa- 
rameters is completely restored after annealing at 400 °C 
whereas the c- lattice parameter appears to be saturating 
at a value near 0.6% larger than that of the undamaged 
sample. 

The superconducting transition temperatures were de- 
termined using a 1% screening criteria in zero field cooled 
DC magnetization measurements (Figure 3a) and an on- 
set criteria in resistivity versus temperature measure- 
ments (Figure 3b). The magnetic transitions for the en- 
tire set are fairly sharp, typically achieving 95% of the 
saturated value 2.5 K below the 1% criteria used for Tc. 
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This suggests a relatively homogeneous damaging of the 
sample. We were unable to obtain reliable transport mea- 
surements on as damaged as well as the 100 °C annealed 
samples due to the 120 °C temperature required to cure 
the silver epoxy contacts. The 30 minute 120 °C cure 
of the epoxy, as well as the possible exotherm associated 
with the epoxy curing, placed the sample into a poorly 
defined annealed state. As shown in figures 3a and 3b, 
superconductivity is restored by the annealing process. 
ATc, defined as the difference between the undamaged 
Tc and that of the annealed sample, monotonically ap- 
proaches zero as the annealing temperature is increased 
(Figure 4). 

The zero field temperature dependent resistivity data 
are plotted in figure 3b. The sample annealed at 150°C 
has a low temperature normal state resistivity slightly 
above 6 cm, which is an order of magnitude larger 
than the undamaged value of 0.2 /if2 cm. After anneal- 
ing at a temperature of 200°C, p increases by another 
four fold. Further increases in the annealing tempera- 
ture result in a monotonic decrease of the normal state 
resistivity, to 1.6 cm at an annealing temperature 
of 500°C. Although we were able to perform a measure- 
ment on only a single sample, this increase in p at 200°C 
is thought to be a real effect, as it also manifests itself 
in normalized p/p{300 K) plots, indicating it is not as- 
sociated with a pathologic geometric problem such as 
cracks. Additionally, un-damaged samples consistently 
showed Tc values near 39 K and normal state resistivi- 
ties in the 0.2-0.4 fj,fl cm range, values in agreement with 
previously published resultsSiS^, suggesting the effect is 
the result of the irradiation and post-exposure annealing 
process. Such an anomalous increase in resistivity as a 
function of annealing temperature was also observed in 
the case of neutron irradiation on carbon doped MgB2 
wire segments^. The initial annealing may result in a 
change in the defect structure. For example, vacancies 
don't necessarily immediately recombine with intersti- 
tials to eliminate defects. Vacancies can initially cluster 
together to form higher order defects, which can be rel- 
atively stablc^^. Thus the initial increase in resistivity 
may result from a reorganization of defects into higher 
order clusters which may enhance scattering. 

Insight into the nature of the defects can be obtained 
by tracking changes in the parameter Ap=p(300 K)-p(40 
K). Figure 3d plots Ap as a function of the transition 
temperature. Ap decreases monotonically as a function 
of the transition temperature. This result is in contrast to 
the results of He ion irradiation on MgB2 thin films^^ , but 
in agreement with the predictions of Mazin et al.^^, who 
claimed that as a material is driven into the dirty limit, 
the conductivity in the pi band decreases, causing normal 
state electric transport to switch from the pi to the sigma 
band. Indeed, for the extreme clean and extreme dirty 
limits they obtain values ranging from near 10 to near 
35 pQ cm respectively, consistent with our experimental 
observations of 6 and 35 pfl cm. 

An alternate framework for understanding this mono- 



tonic change in Ap is to recall that the sample is poly- 
grain with random orientation of the crystallites. Given 
that the electrical resistivity of MgB2 is anisotropic to 
start with, combined with the fact that damage appears 
to more aggressively affect the c-axis periodicity, it is 
possible that for lower Tc values an increasing number of 
grain orientations will present high resitivity paths to the 
current. This will effectively reduce the cross sectional 
area of the wire sample. This would result in a similar 
behavior to what has been found for samples that have 
either poor connectivity or gross impurities on the grain 
boundaries^S. 

The upper critical field was determined using an on- 
set criteria in resistivity versus temperature in applied 
fields up to 14 T (a representative set is shown in fig- 
ure 5a) and, in the case of the 500°C anneal, resistance 
versus field sweeps up to 32.5 T (Figure 5b). The up- 
per critical field curves nest, forming a sort of Russian 
doll pattern, with IIc2(T=0) approximately scaling with 
Tc (Figure 6). The curves for samples annealed at tem- 
peratures up to 200°C do not show any positive curva- 
ture near Tc and are qualitatively similar to single gap 
superconductors with Werthamer, Helfand, and Hohen- 
berg (WIIII)Si like behavior. Experimentally determined 
Hc2(T=0) values for the 150°C, 200°C, and 300°C an- 
neals are 2.9 T, 4.7 T, and 7.3 T respectively. Using the 
formula IIc2(T=0)=0.69TcdIIc2/dT we obtain estimates 
of 2.9 T, 4.3 T and 5.9 T. Thus whereas only in the 
cases of the 150°C and 200°C anneals we can fit Hc2(T) 
with WHH behavior, the deviations increase with the 
annealing temperature, suggesting that the bands may 
become fully mixed only when Tc is suppressed to near 
10 K. Single gap behavior has been inferred from specific 
heat measurements on irradiated samples containing iso- 
topically enriched ^^B which had Tc near 11 K-^. The 
300°C, 400°C, and 500°C anneals exhibit positive cur- 
vature near Tc that is similar to what is found in pure 
MgB2. The 500°C anneal data show that either the un- 
damaged IIc2 (T) is restored or that there a slight increase 
in IIc2(T=0), rising from approximately 16 T in the un- 
damaged case to near 18 T. 

This behavior in Hc2 differs from other types of neu- 
tron damaging studies. For the case of irradiating iso- 
topically enriched Mg^^B2. Putti et. al found a fairly 
substantial increase in Hc2 values when Tc was in the 36- 
38 K ranged. For a sample with Tc = 36.1 K, they re- 
port Hc2(T=12K) of 20.3 T. For unshielded irradiation of 
MgB2 containing natural boron, Eisterer et. al reported 
that, relative to an undamaged sample with a Tc just 
below 38 K, a sample with a suppressed Tc of approxi- 
mately 36 K had a near doubling of the slope dllc2 /dT 
in the linear regime above 2 T^^. In both of these cases 
the fluences used were an order of magnitude lower than 
in our experiment and no post exposure annealing was 
performed. 

Figure 7 presents the critical current densities for the 
series of 24 hour anneals at 5 K and 20 K as determined 
by the Bean critical state model'^'* from magnetic hys- 
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teresis loops. At both 5 K and 20 K, Jc is suppressed for 
samples annealed at 300 °C and below, presumably due 
to the larger reduced temperature (T/Tc) of the mea- 
surements. In field Jc values become enhanced at both 5 
K and 20 K for the 400 °C and 500 °C anneals, and show 
a slight fishtail or second peak behavior. This is believed 
to be the first time such behavior has been observed in 
polycrystalline samples. The fish tail effect has been seen 
in fast neutron irradiation of MgB2 single crystals^!. The 
authors report a second peak in the magnetization hys- 
teresis loops for H_Lab, which became more pronounced 
with increasing fluence up to AV x 10^^ cm~^. No second 
peak was observed for H||ab. Low level thermal neutron 
fluence, of order 10^^-10^^ cm^^, of polycrystalline MgB2 
showed an enhancement in flux pinning, but no second 
peak or fishtail behavior— . The infield enhancement of 
Jc for the sample annealed at 500°C pushes the point at 
which Jc crosses lO'' A/cm^ at 20 K out to near 1.5 T. 
This increase in field roughly doubles the crossing point 
of the undamaged wire, but falls well below best reports 
in literature of 10^ A/cm^ at near 5 T reported for 10 
weight percent addition of SiC to MgB-^. 

Other authors have analyzed critical current densities 
in neutron irradiated samples using a percolative model 
that considers, as free parameters, the upper critical field, 
the absolute value of the critical current density, the 
anisotropy of the upper critical field, and the percola- 
tion thresholcSiiSS*^. In reference^ the authors extend 
the model to include different contributions from grain 
boundary and precipitate pinning. The model yields ex- 
tremely good fits to the data, from which the authors 
conclude that the anisotropy of Hc2 heavily influences 
the fleld in which the material can carry large amounts 
of current''"'^''^^ and defects introduced by the irradiation 
process are responsible for the enhancement of the pin- 
ning strengtlr^^. It is worth noting that these reports do 
not show the significant change in Jc(0) that we observe 
(Figure 7). Presumably two effects are causing the in- 
crease as a function of annealing temperature and hence 
Tc. In using the percolation model, each of the authors 
kept the percolation threshold fixed in the 0.2-0.3 range. 
If we are observing large changes in the connectivity in 
our samples as a function of annealing temperature, as 
inferred from the Ap data, then the initial decrease in 
Jc(0), as observed in the low temperature anneals, is the 
result of the formation of weak links which limit the over- 
all critical current density. As the annealing temperature 
is increased the connectivity between grains is improved, 
thereby returning Jc(0) towards the undamaged value. 
In the case of the 500°C anneal, Jc(0) exceeds that of 
the undamaged sample, suggesting that the defects which 
remain following post-exposure annealing act as effective 
pinning sites. If defects caused by irradiation are 0.1-10 
nm in size, as reported in reference^, then following the 
annealing process, many of the defects should be of or- 
der the coherence length-^, making them effective flux 
pinning centers. 

A plot of Jc as a function of reduced fleld (II/IIc2) 



shows that for the higher temperature anneals, the frac- 
tion of IIc2 in which samples can maintain in excess of 10'^ 
A/cm^ is greatly enhanced (Figure 7c and 7d). Such en- 
hancement is consistent with an increase in the strength 
of defect pinning'^®. That the low temperature anneals 
all show Jc dropping below 10^ A/cm^ at similar reduced 
fleld values suggests the defects are ineffective in enhance 
critical current densities in highly disordered samples. 
This could be the result of, for example, overlap of de- 
fects, making them to large to be effective flux pinning 
sites. 



C. Variable Time Anneals on Samples Exposed to 
a Fluence of 4.75 x 10^** cm"^ 

The annealing time for the set of 4.75 x 10^^ cm~^ 
fluence samples was varied to further probe the charac- 
teristics of damage induced by neutron irradiation. An 
extensive set of anneals was carried out at 300°C, consist- 
ing of 1/3, 1, 3, 6, 24, 96, and 1000 hours. For annealing 
temperatures of 200°C and 400°C, measurements were 
done on 1, 24, and 96 hour anneals. Two anneals were 
performed at 500°C, one for 24 hours and the other for 
1000 hours. 

The magnetic transitions for the entire set of samples 
annealed at 300°C are shown in flgure 8. Annealing at 
this temperature for only 0.33 hours raised Tc from below 
5 K to slightly above 19 K. Therefore, the defects causing 
the suppression of superconductivity must have a fairly 
low activation energy. For defects which can be annealed 
by a single activated process with a constant activation 
energy, the rate of change of the defect concentration is 
given bjii^: 

dn/dt -F(n)Xoe"-^°/''^^ (2) 

where n is the defect concentration, F{n) is some con- 
tinuous function of n, Kq is a constant, Ea is the ac- 
tivation energy, fc^ is Boltzman's constant, and T is 
the temperature. If we assume random diffusion, then 
F{n) = n, and the defect concentration decreases expo- 
nentially with time: 

n = noe~^* (3) 

where no is the initial defect concentration, C — 
Kq g-^i/'^sT ^ j.g^-(-g constant, and t is time. 

While po can be taken as one measure of the defect 
concentration, we have already shown non-nionotonic be- 
havior in Po as a result of the annealing process. An- 
other variable that can be used to defect concentration 
is Tc, which systematically varies with both annealing 
temperature and annealing time. We therefore take ATc 
as a measure of the defect concentration, although it 
should be reiterated that the nature of the defects may 
be changing as we increase the annealing temperature 
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and time. Such linearity is not uncommon and can be 
related to Abrikosov-Gorkov-like theory of pair break- 
ing in anisotropic superconductors 4^ It should be noted, 
though, that we do this simply in an attempt to extract 
an approximate value of the activation energy. 

We see a exponential decay behavior in ATc for sam- 
ples annealed at 200°C, 300°C and 400°C (Figure 9). 
Since we only have two samples annealed at 500°C, which 
are longer time anneals and thus presumably would fall 
further out on an exponential tail, we can not extract a 
meaningful value from an exponential fit of these data. 
For samples which do show an exponential behavior in 
the decrease of ATc (and by presumption) the defect 
density as a function of annealing time the activation 
energy, assuming a single activation process, can be de- 
termined by the so-called cross-cut procedure^ii. This 
involves comparing the annealing time for which differ- 
ent temperature anneals reached the same defect density, 
i.e. ATc- If a ATc is reached by annealing at a tempera- 
ture Ti for a time ti and by annealing at a temperature 
T2 for a time t2, then the activation energy is related to 
these quantities by: 

ti Ea , 1 1 , , , 

Unfortunately, no overlap region in Tc exists for our 
set of data. To achieve an overlap region in the exper- 
imental data we'd need to shorten the annealing times 
below 0.33 hours or extend them to well beyond 1000 
hours. Short time anneals, less than 0.33 hours, are not 
feasible as we'd be unable to ensure thermal equilibrium 
of the samples within the furnace for such a short time. 
Extending the annealing an additional order of magni- 
tude, to IC* hours, is simply not practical. We can still 
estimate the activation energy by extrapolating the ATc 
curves at each temperature to shorter and longer times so 
as to create overlap regions (Figure 9). Calculating Ea 
from the overlap between the 200°C and 300°C curves 
and between the 300°C and 400°C curves yields values 
of 1.90 eV and 2.15 eV respectively. 

An alternative approach to determining the activation 
energy is to use the ratio of slopes method"^^. Exper- 
imentally a set of samples is annealed isothermally for 
different times at a temperature Ti and an identical set 
of samples is annealed at a different temperature T2. The 
differing temperatures result in different time evolution 
of ATc, and hence, different slopes, dATc/dt. The ratio 
of the slopes for the two different temperatures at the 
point where both annealing temperatures have yielded 
the same ATc, is related to the activation energy by: 

dti dt2 Kb J-i J-2 

where the subscripts 1 and 2 refer to temperatures Ti and 
T2. By comparing the slopes from the 200''C and 300°C 
anneals and those from the 300°C and 400°C anneals, we 



obtain estimates of 1.07 eV and 1.63 eV, respectively. It 
should be noted that, while these estimates were made 
using real data, rather than extrapolations as in the pre- 
vious calculation, there is inherent inaccuracy in such 
a calculation due to the lack of a true overlap in ATc. 
Additionally, the low density of data points limits our 
ability to accurately determine the linear slope from the 
ATc versus time curves where ATc tends to decay expo- 
nentially. 

Depending on the calculation method and data set 
used, we obtain activation energies ranging from 1.07 eV 
to 2.15 eV. It is likely that some of the variation is real, as, 
in these heavily damaged samples, there exist both point 
defects and defect complexes. The annealing of point de- 
fects is expected to have a lower activation energy than 
the dissolving of the defect complexes. While we can not 
assign definitive values for the activation energies of these 
two processes, merely stating a single activation energy 
hides some of the rich complexity underlying the anneal- 
ing process in these heavily damaged samples. It should 
be noted that the activation energies in these neutron ir- 
radiated MgB2 samples are the same order of magnitude 
as those for annealing quenched in defects out of gold^S. 
That our samples yielded a relatively small spread in ac- 
tivation energies and were comparable in magnitude to 
values associated with the annealing of simple defects 
in other metals suggests the defects within the neutron 
irradiated MgB2 samples are being annealed by single 
activation processes. 

For the set of wires annealed at 300°C, critical cur- 
rent densities at T = 5 K and low fields approximately 
scale with annealing time and hence Tc (Figure 10). The 
field at which Jc drops below 10'' A/cm^ increases by 
approximately a factor of two when the annealing time 
is increased from 1 hour to 24 hours. Extending the an- 
nealing time further to 1000 hours has little effect on low 
field Jc values. For annealing times of 6 hours and shorter 
(for the sake of clarity, data for 1/3, 3, and 6 hours are 
not shown) Jc was found to monotonically decrease as a 
function of applied field. After a 24 hour anneal, Jc be- 
gins to flatten above 3 T at a value of approximately 500 
A/cm^. Extending the annealing time to 96 and 1000 
hours results in the emergence of a clear second peak at 
an applied field slightly above 3 T. The curves for the 96 
and 1000 hour anneal are virtually identical and reflect 
the small 1.5 K or 6% increase in Tc resulting from the 
order of magnitude longer annealing time. 



D. Variable Fluence Levels, 24 Hour Anneal 

We annealed samples from each of the damaged levels 
for 24 hours at 200°C, 300°C, 400°C, and 500°C. The 
(002) and (110) x-ray peaks for samples of all four dam- 
age levels annealed at 300°C for 24 hours are given in fig- 
ure 11. Annealing at 300°C for 24 hours did not restore 
long range order along the c-axis in the samples exposed 
to the two highest fluence levels of 1.43 x lO"'^ cm~^ and 
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1.90 X 10^^ cin~^. For these 24 hour anneals, up to an 
anneahng temperature of 400°C the (002) peak main- 
tains a FWHM above 1°, which corresponds to a struc- 
tural coherence length of approximately 1000 A (Figure 
12). Using Hc2(T=0) to determine the supercoducting 
coherence length, one obtains for pure MgB2 a super- 
conducting coherence length near 50 A^^ . Measurements 
on the coherence length within the tt band only yield a 
larger value of approximately 500 A^. While there is 
shorter range structural order along the c- direction in 
these heavily irradiated samples, the structural coher- 
ence length is still at least a factor of two larger than 
the superconducting coherence length, regardless of the 
estimate of the superconducting coherence length used. 
It is therefore not unexpected that superconductivity ex- 
ists in samples with degraded long range order. Only 
after the temperature reaches 500°C does the correlation 
length along the c- direction exceed 1000 A (Figure 12). 
The relative shift of the a- and c- lattice parameters for 
the various heat treatments is plotted in figure 13. 

Magnetization measurements were performed on all of 
the annealed samples (Figure 14a). It should be noted 
that for the fluences greater than 4.75 x 10"'^^ cm~^, as 
damaged samples showed no signs of superconductivity 
down to 2 K. A summary of the transition temperatures 
for each of the annealed samples, given in terms of ATc, is 
plotted in figure 14b. In general terms, the higher fluence 
levels lead to lower superconducting transition tempera- 
tures for a given temperature post-exposure anneal. For 
a given fluence level, higher annealing temperatures yield 
higher Tc values. That is, the samples exposed to higher 
fluence levels behave in a qualitatively similar manner to 
lowest level, but the added exposure leads to increased 
defect densities which manifest themselves in terms of 
lower transition temperatures for a given annealing pro- 
file. 



E. Long Time Annealing Studies 

The conversion of ^"^B to ^Li through neutron absorp- 
tion and subsequent alpha decay introduces the possibil- 
ity of observing the effects of lithium doping MgB2. In 
order to distinguish the effects of Li doping from those 
associated with structural defects introduced through in- 
elastic collisions between neutrons and the underlying 
lattice it is necessary to minimize the density of these 
defects. As shown, by annealing for long times at high 
temperatures the number of structural defects can be sys- 
tematically reduced. Since Li can't be annealed away, the 
resultant superconducting and normal state properties 
should increasingly reflect the effects of Li doping. The 
density of Li atoms produced through the transmutation 
of boron can be estimated from the formulcii^: 

riLi = riBsf (6) 
where n^i is the density of Li atoms, is the density 



of B atoms, s is the absorption cross section, and / is 
the fluence level. Computing the corresponding atomic 
percentages yields an increase from 0.37% to 1.48% as the 
fluence in increased from 4.75 x 10^^ cm^^ to 1.90 x 10^^ 
cm" 2 (Table 1). 

We therefore annealed all four of the damage levels for 
1000 hours at 500°C. Normalized magnetization curves 
for this series are plotted in figure 15a. In each case the 
transition temperatures have increased relative to sam- 
ples annealed at 500 °C for 24 hours as can be seen by 
plotting ATc versus time (Figure 15b). Using a linear 
fit on the semi- log plot, the resultant slopes are much 
lower than was seen for the 4.75 x 10^^ cm~^ fluence level 
samples annealed for various times at 200°C, 300°C, and 
400°C (Figure 9), indicating ATc is beginning to satu- 
rate above 24 hours. Since we don't have any interme- 
diate time points, we can't determine if we've achieved 
fully saturated ATc values at 1000 hours, but we can take 
the 1000 hour anneals at 500°C as an upper limit on the 
effects of lithium doping. 

Transport measurements were performed in order to 
determine the normal state resistivity and temperature 
dependance of Hc2 for these samples. We were unable to 
contact the samples exposed to the two highest dose lev- 
els. The normal state resistivity (T=40 K) for the sample 
exposed to a fluence of 4.75 x 10^* cm~^ and annealed for 
1000 hours was approximately 5.3 fifl cm, which is more 
than three times the 1.6 cm measured on the sample 
annealed for 24 hours. The 1000 hour anneal sample also 
had a lower residual resistivity ratio (RRR), 3.3 versus 
5.9, indicating the increase in po isn't merely an artifact 
associated with possible geometric effects (cracks within 
the sample). The 9.50 x 10^* cm~^ fluence level sam- 
ple, which was annealed for 1000 hours, had a normal 
state resistivity of 6.6 ^fl cm and RRR=5.0. Karkin et 
al. saw increases in p for annealing temperatures above 
300°C, which they attributed to changes in inter-grain 
transport^^. It is possible that by annealing at an ele- 
vated temperature for 1000 hours we have degraded the 
inter-grain connectivity in some fashion, leading to the 
observed increases in resistivity. 

The temperature dependance of Hc2 is plotted in flg- 
ure 16. It was found that for the 4.75 x 10^* cm~^ flu- 
ence level i?c2(T) data for a 1000 hour anneal sample 
was comparable to that of the 24 hour anneal sample, 
both showing possible, slight enhancement relative to the 
pure sample. In the case of the 1000 hour anneal on the 
9.50 X 10^^ cm~^ fluence level, Hc2 had some spread in 
the data, but extrapolated to 18-19 T at zero Kelvin. 
The upper critical flelds of the two lowest fluence levels 
annealed at 500°C for 1000 hours appear to be similar, 
with slight differences arising due to inherent sample to 
sample variation. 
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IV. DISCUSSION 

The initial irradiation of MgB2 wire segments results 
in an increase in the size of the unit cell and suppres- 
sion of the superconducting transition temperature. Post 
exposure annealing tends to return both the lattice pa- 
rameters and Tc towards their undamaged values. If 
the superconducting properties of the neutron irradiated 
samples were purely a result of changes in the unit cell 
dimensions, correlations should exist between the super- 
conducting properties of neutron irradiated MgB2 and 
pure MgB2 placed under external pressure. Since an ex- 
pansion is qualitatively analogous to an effective negative 
pressure one would anticipate the the changes in T^, to 
be an extension to negative pressure of the results at- 
tained for the application of positive pressure. Applica- 
tion of external pressure has been shown to compress the 
unit cell and suppress T,~^. Thin films grown epitaxially 
on (0001) sapphire substrates exhibited an enhanced Tc 
above 41 K which was attributed to tensile strainp^. In 
the case of neutron irradiated MgB2 , the expansion of the 
lattice parameters coincides with a decrease in Tc and the 
evolution of Tc as a function of Aa, Ac, and V/Vq be- 
haves differently than MgB2 under pressure (Figure 17). 
In the case of externally applied pressure, it is believed 
that the changes in the frequency of the F,2g phonon mode 
are responsible for the suppression of Ttp2i^. With such 
dramatically different behavior between the neutron irra- 
diation and pressure results, it is clear that the changes 
in Tc can not simply be linked to changes in the unit cell 
volume. 

In addition to the structural changes, neutron irradia- 
tion also introduces a chemical impurity into the system 
through the ^°B neutron capture and subsequent alpha 
decay to ^Li. The amount of lithium produced in this 
manner is, however, quite small and for all but the 1000 
hour anneals at 500°C, we can not sort the possible effects 
of such from the effects of the structural perturbations. 
The atomic percentage of B converted to Li through nu- 
clear processes for all four exposure levels is estimated 
to be on the order of 1% (Table 1). Therefore, for lower 
temperature and/or shorter time anneals the effects of 
Li production are presumably insignificant next to the 
changes resulting from structural damage caused by in- 
elastic collisions between the fast neutrons, emitted alpha 
particles, and recoiled ^Li atoms and the underlying lat- 
tice. It is worth noting, though, that whereas low level Li 
substitution for Mg was proposed as a possible route to 
increasing Tc^-, our data show suppressed Tc values for 
our 1000 hour annealed samples. This does not preclude 
the possibility that Li substitution could, under other cir- 
cumstances, raise Tc, but it certainly does not support 
it. It should be noted though that this form of lithium 
doping is far from the ideal one. In the case of trans- 
mutation of ^°B the resulting sample is MgLij;B2-2: with 
boron vacancies and no clear site for the lithium. The 
more desirable form of lithium doping is Mgi_a;Lia:B2 
with the lithium substituting for the magnesium and no 



disruption of the boron sublattice. 

The samples annealed at 500°C for 1000 hours suggest 
that at low levels, Li doping has little or no effect on 
Hc2- The 9.50 x 10^^ cm~^ fluence sample has twice the 
amount of Li of the 4.75 x 10^^ cm~^ fluence sample, yet 
their Hc2 (T=0) values are approximately equal and only 
slightly differ from the un-damaged case. If the 18 T IIc2 
value seen for the 4.75 x 10^* cm^^ fluence level is truly 
a 2 T enhancement relative to the un-damaged sample 
and is the result of Li doping, then we would expect 
the 9.50 X 10^^ cm'^ fluence level to exhibit a Hc2{T = 
0) near 20 T, which we did not observe. Since we can 
not insure that we have fully annealed out all of defects, 
it is possible that this slight increase in HcziT = 0) is 
a result of scattering associated with structural defects. 
In this case Hc2{T = 0) should be determined more by 
Tc than by the particular fluence level or post exposure 
annealing profile. The Tc values for the 4.75 x 10^^ cm~^ 
fluence level annealed at 500°C for 24 and 1000 hours 
and that of the 9.50 x 10^* cm~^ fluence level annealed 
at 500°C for 1000 hours are all within half a degree of 
one another. This suggests that Hc2{T = 0) is controlled 
more by scattering associated with residual defects rather 
than any inadvertent Li doping. Since these low levels 
of Li do not appear to have any major impact on the 
superconducting properties, we can limit the discussion 
to possible influences of disorder, scattering, and possible 
changes in the Fermi surface. 

The evolution of the superconducting transition tem- 
perature as a function of the unit cell dimensions (Aa, 
Ac, and V/Vq) shows definite trends and is in good 
agreement with the results of other neutron irradiation 
studiesi^iS^ (Figure 18). Although there is considerable 
spread in the data, all three reports show that Tc tends 
to decrease with both Aa and Ac and hence V/Vq. These 
data show that, for neutron damaged samples, there is 
some correlation between the unit cell dimensions and 
the superconducting properties. They do not, however, 
uniquely determine if the changes are a result of changes 
in the Fermi surface, perhaps due to a repositioning of the 
atoms, or if they are due to an introduction of additional 
scattering centers. In the case of substantial neutron ir- 
radiation induced damage, where Tc is below 10 K, NMR 
measurements indicate the suppression of superconduc- 
tivity is the result of a decrease in the density of states of 
the boron px,y orbitala^Si. This technique has been suc- 
cessful in experimentally determining a decrease in the 
density of states in Mgi_2;AUB2 and A\B2^^. It should 
be noted however, that measurements of the nuclear spin 
relaxation rate, Tj~^, on ^^B yielded comparable TiT val- 
ues, which are directly related to the density of states at 
the Fermi surface, for both a neutron irradiated sample 
which had a Tc near 7 K"*^ and a sample with 30% alu- 
minum substituted for magnesium which had a Tc near 
25 K*'^ . That the neutron damaged sample exhibits a 
dramatically lower Tc despite having virtually the same 
density of states at the Fermi surface suggests additional 
mechanisms act to suppress superconductivity in neutron 
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irradiated MgB2. This notion is supported by the fact 
that for the 24 hour anneals at 300°C, 400°C, and 500°C, 
the relative changes in the lattice parameters are fairly 
small. This should only contribute to minor changes in 
the Fermi surface, yet Tc is suppressed to near 25 K in 
the case of the 300°C anneal. 

Although there exist correlations between our data and 
literature reports on the evolution of Tc and the lattice 
parameters, the evolution of i?c2 [T = 0) varies depend- 
ing upon the irradiation conditions. This present work 
on heavily irradiated MgB2 containing natural boron fol- 
lowed by post exposure annealing led to little or no en- 
hancement of the upper critical field. In contrast others 
have reported enhancements of Hc2 for low fluence levels 
on MgB2 containing either natural boron or isotopically 
enriched ngis ^o ^^qW fg^g^ neutron irradition of nat- 
ural boron containing MgB2'^"'^. 

In order to better understand this conspicuous differ- 
ence it is useful to review what happened during dam- 
age and subsequent annealing. When a neutron capture 
and subsequent alpha decay event occurs, the resultant 
damage is a large cluster of dislocations. Primary knock- 
on events cause a cascade of displacements which can 
be spread over a distance as large as a hundred atomic 
distancesii. 

As mentioned previously, Tc suppression in as dam- 
aged samples containing natural boron is much more 
rapid than those with either isotopically enriched ^^B or 
those which have been shielded from low energy neutrons. 
Therefore it is presumably the clusters associated with 
the neutron capture event which are primarily responsi- 
ble for the suppression of Tc. These clusters presumably 
have low activation energy and can be largely repaired by 
the annealing process, which explains why we saw such a 
rapid increase in Tc for short time and low temperature 
anneals. In order to see an enhancement in Hc2 we need 
to anneal for a sufficient time and temperature so as to 
increase Tc to a level where we aren't being limited by 
a low transition temperature. Since the defect clusters 
from the alpha decay of ^^B to ^Li are presumably large, 
and we are relying on random diffusion to recombine va- 
cancies and interstitials, higher reaction temperatures are 
necessary to restore Tc- 

By going to higher annealing temperatures we begin to 
repair defects with higher activation energies. Therefore, 
regardless of the actual energy associated with repair- 
ing the defects responsible for the enhancement of Hc2, 
we are most likely annealing many of them away while 
restoring Tc- 

If it is these defects which are responsible for Hc2 en- 
hancement then it is the annealing process which causes 
the different development of Hc2 we observe. 

Based on this analysis it seems likely that a study of 
wire samples that have been exposed to much smaller 
fluences of neutrons (e.g. an exposure level that yields 
as damaged samples with a transition temperature in the 
range 30 < Tc < 37) should yield different, and higher, 
Hc2iT) data. 



The transmutation of ""^"^B following the absorption of 
a neutron results in the emission of an alpha particle of 
energy 1.7 MeV, which is quite comparable in energy to 
the 2 MeV He ion irradiation of MgB2 films done by 
Gandikota et al.^*'^^. One might therefore expect rather 
similar results between the two studies. Comparison of 
the results indicates these two different types of irradia- 
tion give rise to different superconducting and normal 
state properties. Gandikota and coworkers irradiated 
with heavy ion fiuences up to 1.2 x 10^^ cm~^, which was 
sufficient to suppress Tc to below 5 K^^. Unlike the case 
of heavy neutron irradiation presented here, He"*""*" irradi- 
ation led to a slight enhancement in Hc2 [T = 0) for sam- 
ples exposed to lower fluence levels, where Tc was only 
suppressed by a few degrees^, analogous to the results 
obtained for lower fiuence level neutron irradiation^^'^Si. 
Hc2{T = 0) did tend to scale with Tc for samples with 
Tc below 30 K and, as with neutron irradiated samples, 
Hc2 {T) begins to exhibit more WHH like behavior when 
Tc in the vicinity of 10 K. Additionally, Tc could read- 
ily be restored by annealing at temperatures as low as 
100°C''^. However, x-ray measurements indicated that 
the He ion irradiation did not change the c-lattice pa- 
rameter, even for the highest fluence level2&. Thus, He"'"'" 
irradiation presumably creates more point defects rather 
than large defect clusters which are likely to arise from 
neutron irradiation. He ion irradiation therefore presents 
yet another different avenue to systematically change the 
scattering levels in MgB2. 

Damage induced by neutron irradiation is fundamen- 
tally different from doping with aluminum or carbon. 
Both carbon and aluminum enter the structure and are 
believed to act as point defects in addition to electron 
doping the system. Arguments based on the tempera- 
ture dependence of the anisotropy of the upper critical 
field suggest both carbon and aluminum doping increase 
scattering within the tt band relative to the a band, with 
the effect is much more pronounced in the case of carbon 
dopingi. It should be noted that other researchers have 
concluded that changes in the temperature and magnetic 
field dependence of the thermal conductivity, k(T,H), for 
carbon substituted single crystals are consistent with car- 
bon doping resulting in an enhancement in intra-cr-band 
scattering^'^. Although there exists some debate as to the 
exact nature of the enhancement in scattering in carbon 
doped MgB2, it is clear that the development of Hc2 is 
dominated by scattering effects for carbon substitutions 
and and Fermi surface changes effects for aluminum sub- 
stitution. For neutron irradiated samples, as discussed 
above, the effects of the introduction of Li through nu- 
clear processes are negligible and the changes in the su- 
perconducting properties are due to an increase in inter- 
band scattering and a decrease in the density of states. 
Insights can be made into the effects of neutron irradiated 
samples by comparison to carbon and aluminum doping. 

The different development of interdependencies of 
jjll ab^ Tc, and po in aluminum doped^, carbon dopedSS*^, 
and neutron irradiated samples, those exposed to a flu- 
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ence of 4.75 x 10^* cm~^ and annealed for 24 hours, are 
plotted in figure 19. Figure 19a shows the evolution of 
the transition temperature as a function of residual re- 
sistivity for these three types of perturbations. The sup- 
pression of Tc with po for low resistivity values is quite 
similar for the neutron damaged and aluminum doped 
samples; drops rapidly with increased scattering. In 
contrast, the suppression of Tc as a function of resistivity 
in carbon doped samples is quite gradual. This is consis- 
tent with both aluminum doping and neutron damaged 
samples having more interband scattering than carbon 
doped compounds. 

Figure 19b plots Hc2{T = 0) as a function of po. Here 
all three perturbations behave uniquely. The aluminum 
and neutron damage samples show a decrease in Hc2, 
whereas the carbon doped samples show a dramatic in- 
creases. The decrease in Hc2 is more rapid in the alu- 
minum doped samples than in neutron damaged samples. 
This can also be seen by plotting Hc2 as a function of Tc 
(Figure 19c). 

Direct comparison between the evolution of Tc and 
Hc2{T — 0) for carbon doping and neutron irradiation 
shows the scattering associated with each type of per- 
turbation is different. Within the model proposed by 
Gurevichii, enhancements in Hc2 result from differences 
in the relative strength of scattering within each band, 
whereas the suppression of Tc is a function of scatter- 
ing between the bands. The neutron irradiated samples 
on these plots with a Tc greater than 25 K are samples 
which were annealed at 300°C, 40Q°C, and 500°C. In 
each of these samples, annealing has reduced the a-lattice 
parameter to nearly the undamaged value (Figure 2b). 
With the defect structure lying between nearly undam- 
aged boron planes, one would not expect the scattering 
within the a band to be substantially affected. If the scat- 
tering was confined primarily to the 3-D tt band, the re- 
sultant differences in intraband diffusivity values should 
manifest themselves in terms of enhanced Hc2- Since no 
significant enhancement is seen, we believe the scatter- 
ing is primarily interband scattering and contributes to 
the suppression of Tc- The range in Tc for V/Vq ^ 1 
in post exposure annealed samples (Figure 18c) supports 
the notion that neutron irradiation increases interband 



scattering. 

The suppression in Tc in neutron irradiated samples 
is the result of a combination of a decrease in the den- 
sity of states at the Fermi surface and an increase in 
interband scattering. In the case of aluminum doping, 
the suppression of Tc is believed to be primarily a re- 
sult of a changes in the Fermi surface^. Samples of neu- 
tron damaged MgB2 exhibit higher Hc2 values than alu- 
minum doped MgB2 samples with similar Tc values (Fig- 
ure 19c). This suggests that, for neutron irradiated sam- 
ples, changes in the density of states suppress Hc2 to a 
greater degree than interband scattering. 

V. CONCLUSIONS 

We systematically studied the effects of neutron fluence 
level, annealing temperature, and annealing time on the 
superconducting and normal state properties of MgB2. 
As damaged samples showed an anisotropic expansion of 
the unit cell and a suppression of Tc to below 5 K. Defects 
introduced by the irradiation process had relatively low 
activation energy and hence much of the damage could be 
repaired by post-exposure annealing. Hc2 values tended 
to scale with Tc and the evolution of Hc2 with po was 
dramatically different than in the case of carbon doping. 
Lithium produced through nuclear processes appeared to 
have no significant effects on the superconducting prop- 
erties. We attribute the changes in the superconducting 
properties of neutron irradiated samples primarily to an 
increase in scattering between the two bands and possible 
changes in the density of state at the Fermi surface. 
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FIG. 1: (a) (002) and (110) X-ray peaks for un-damaged, 
4.75 X 10^^ 9.50 X 10^^ 1.43 x 10"\ and 1.90 x 10^'' mr^ 
flucncc as-damaged samples. The highest two fiuenco levels 
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FIG. 3: (a) Normalized magnetization and (b) resistivity 
curves for the set of 24 hour anneals on samples exposed to a 
fluence of 4.75 x 10^* cm~^. (c) The normal state resistivity 
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nealing temperature of 200 °C, then decreases approximately 
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FIG. 9: (Color online) Semi-log plot of the change of the 
superconducting transition temperature as a function of an- 
nealing time at various annealing temperatures. All samples 
shown were exposed to a fluence level of 4.75 x 10^* cm~^. Us- 
ing a linear fit over up to three decades yields non-systematic 
values in the rate constants for the 200°C, 300°C, and 400°C 
annealing temperatures. The activation energy is estimated 
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FIG. 11: (002) and (110) X-ray peaks for all four exposure 
levels annealed at 300° C for 24 hours. The highest two levels 
continue to show a substantially broadened (002) peaks indi- 
cating a degradation of long range order along the c-direction. 
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FIG. 12: (Color online) Evolution of the (002) and (110) X- 
ray peaks for the 1.43 x lO^'' and 1.90 x 10^^ cm~^ exposure 
levels annealed for 24 hours with the annealing temperature 
increasing from 300°C to 500°C. In both cases, the (002) at- 
tains a FWHM less than 1° 20 only after the annealing tem- 
perature reaches 500°C. Note: intermediate temperature an- 
neals at 400° C were performed for longer times to show the 
persistence of the broadening for T< 400°C. 
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FIG. 13; (Color online) Calculated lattice parameter shifts 
for all four damage levels annealed at temperatures up to 
500°C. Closed symbols represent Aa, and open symbols are 
Ac. Lines serve as guides to the eye. 
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FIG. 14: (a) Normalized magnetic transitions as a function of 
annealing temperature and fluence level for all four damage 
levels. The time for each anneal was 24 hours. The curves are 
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1.2 units for ease of reading, (b) ATc values, as determined by 
a 1% screening criteria, as a function of annealing temperature 
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exposed to a fluence of 9.50 x 10^* cm~'^ are shown. Both 
damage levels show a possible enhancement relative to the 
undamaged sample, but there is some spread in the data as 
illustrated by the two 9.50 x 10^^ cm~'^ fluence level samples 
shown. 




-0.06 



-0.05 



-0.04 



-0.03 



-0.02 



-0.01 



0.00 



0.01 



Aa (Angstroms) 



40 




-0.10 -0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 0.08 



Ac (Angstroms) 



41 



40 



35 - 



30 - 



25 I- 



15 - 



10 - 



5 - 



★ 


Pure 


• 


Neutrons 


A 


Pressure 



X 



X 



X 



X 



(c) 

— 1— 



X 



X 



X 



0.84 0.86 0.88 0.90 



0.92 0.94 

v/v„ 



-r^ — r 



X 



_i_ 



0.96 0.98 1.00 1.02 



FIG. 17: (Color online) Comparison of the different develop- 
ment of Tc with (a) Aa, (b) Ac, and (c) unit cell volume for 
neutron irradiated MgB2 and MgB2 under external pressure. 
Pressure data are recreated from reference—. 
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FIG. 18: (Color online) Comparison of the different develop- 
ment of Tc with (a) Aa, (b) Ac, and (c) unit cell volume for 
neutron irradiation from different groups. Data includes re- 
sults from references—*—. In parts (a) and (b) the error bars 
on the pure sample represent typical experimental error in de- 
termining lattice parameters. In the part (c), the error bars 
represent the propagation of the error in the lattice param- 
eters to the calculated V/Vo value. All of our data in these 
plots should be considered to have comparable uncertainty. 



45 



45 
40 
35 
30 

o 25 - 
20 - 
15 - 
10 



Pure 



h T 



(a) 







• n-damaged 
^ carbon 
T aluminum 



10 15 20 
cm) 



25 



30 



40 



(b) 



"I ' 1 ' 1 ' 1 ' r 



30 - 



• n-damaged 
^ carbon 
▼ aluminum 



20 - 



Pure 



10 - 







I I I I I I I I I I I 

5 10 15 20 25 30 



p„ (nn cm) 



47 



40 







(c) 







1 ' 


1 ' 


1 ' 1 
▲ 




• 


n-damaged 








A 


carbon 




A 


- 


T 


Aluminum 






— 






▲ 


A 










A 


— 








* * Pure ' 










w 










T 








• 








• 

1 


1 


T 

1 . 1 



10 



20 



30 



40 



T (K) 



FIG. 19: (Color online) Interdependencies of Tc, H^2iT=0), 
and po for carbon doped—*—, aluminum doped single 
crystals-, and neutron irradiated samples, (a) Tc vs. po, 
(b) Hc2(r = 0) vs. po, and (c) H^iiT = 0) vs. T,. 
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TABLE L Estimated fluences for each of the four exposure 
times and calculated atomic percentage of boron transmu- 
tated to lithium. 



